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1.  Abstract 

This  work  combines  two  efforts  that  have  been  on-going  in  the  Air  Force  Research  Laboratory.  The  first  effort  is  the 
work  of  the  Design  Modeling  Group  which  is  using  the  Adaptive  Modeling  Language  (AML)  to  develop  an  object- 
based  technology  assessment  system.  One  of  the  advantages  of  this  effort  is  that  the  analytical  tools  (legacy  codes)  that 
are  used  may  be  integrated  into  the  AML  environment  one  at  a  time,  thereby  providing  many  intermediate  products 
each  of  which  may  enhance  productivity.  So  while  this  effort  is  still  in  its  infancy,  a  graphical  user  interface  to  the 
Automated  STRuctural  Optimization  System  (ASTROS)  has  been  developed  and  is  used  in  this  effort  to  develop  a 
parametrically  defined  finite  element  model.  The  second  effort  is  the  work  on  designing  flexible  wings  without  ailerons 
for  performing  roll  maneuvers.  Here  a  flexible  wing  is  optimized  for  weight  with  constraints  on  strength  for  a  9-g 
symmetric  pull-up  maneuver  at  M=0.85  at  30  psi  and  on  the  modal  frequency  distribution.  Differential  wing  twist  and 
camber  is  achieved  by  providing  a  system  of  actuating  elements  distributed  within  the  internal  substructure  of  the  wing 
to  provide  roll  control  power.  The  modal  approach  was  used  to  develop  equilibrium  equations  for  the  steady  roll 
maneuver  of  a  wing  subjected  to  aerodynamic  loads  and  actuating  forces.  The  distribution  of  actuating  forces  to  achieve 
a  specified  flexible  roll  rate  was  determined  by  using  an  optimal  control  design  approach.  Here,  full-scale  wings  are 
considered  for  the  assessment  of  strain  energy  as  a  measure  of  the  necessary  power  required  to  produce  the 
antisymmetric  twist  and  camber  deformation  to  achieve  the  required  roll  performance.  The  wings  used  here  break  the 
paradigm  of  the  low  aspect  ratio  wings  typical  for  supersonic  vehicles.  The  purpose  of  this  is  to  enhance  the  flexibility 
of  the  wing  and  to  reduce  the  control  power  needed  to  twist  the  wing. 


2.  Introduction 


The  purpose  of  the  technology  assessment  system1  that  is  being  developed  in  the  Adaptive  Modeling  Language  (AML) 
is  to  look  at  how  the  new  technologies  and  concepts  affect  aircraft  design  and  cost.  New  technologies  can  be  assessed 
as  enhancements  to  existing  air  vehicles  or  as  concepts  for  new  vehicles.  Several  areas  of  the  technology  assessment 
system  that  are  currently  under  development  include  a  new  paradigm  for  conceptual  design  tools,2  conceptual  design 
tools  for  hypersonic  vehicles,  a  graphical  user  interface3  for  the  Automated  STRuctural  Optimization  System 
(ASTROS),4  automatic  differentiation5  and  optimization.6  This  work  uses  the  graphical  user  interface  for  ASTROS  to 
develop  an  initial  finite  element  model  for  ASTROS  and  a  second  finite  element  model  based  on  the  first  model  with 
some  minor  planform  changes.  Using  the  parametric  relations  defined  for  the  first  model,  a  second  model  is  generated 
rather  quickly  and  with  minimal  effort. 


Each  of  the  models  generated  for  ASTROS  are 
optimized  for  minimum  weight  and  the 
distribution  of  actuating  forces  to  achieve  a 
specified  flexible  roll  rate  is  determined.  The 
purpose  of  this  portion  of  the  effort  is  to  assess 
the  potential  of  recovering  some  of  the  roll 
effectiveness  that  is  lost  due  to  the  flexibility  of 
the  wing  (see  Fig.  1)  by  using  actuators  to 
control  the  deformation  of  the  wing.7  ASTROS 
uses  the  Unified  Subsonic  and  Supersonic 
Aerodynamic  Analysis  (USSAERO)8  program 
to  compute  the  aerodynamic  loads  on  the 
aircraft  wing.  This  approach  uses  a 
superposition  of  vortex  singularities  applied  to  a 
discrete  number  of  aerodynamic  panels  to 
calculate  the  discrete  pressure  distribution  over 
the  wing  surface. 


Figure  1 .  Effect  of  flexibility  on  the  roll  rate  of  a  typical  wing  for 
various  dynamic  pressures. 
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The  traditional  method  of  using  ailerons  to  produce  a  rolling  moment  is  effective  at  low  dynamic  pressures  where  the 
wing  is  essentially  “rigid.”  However,  at  high  dynamic  pressures,  where  the  deformation  of  the  wing  is  significant,  the 
roll  rate  is  reduced  by  a  detrimental  aerodynamic  twisting  moment  resulting  from  the  trailing  edge  control  surface 
rotation. 

3.  The  Adaptive  Modeling  Language  and  the  Graphical  User  Interface  to  ASTROS 

A  brief  description  of  the  Adaptive  Modeling  Language  (AML)  is  in  order  to  help  understand  some  of  the  benefits  of 
the  Graphical  User  Interface  (GUI)  that  was  developed  for  ASTROS.  AML  is  an  object-oriented,  knowledge-based, 
programming  language  developed  by  TechnoSoft,  Inc.  AML  has  several  unique  capabilities  which  make  it  an  excellent 
language  for  designing  engineering  systems.  These  capabilities  include  dependency  tracking,  demand-driven 
calculations  (lazy  evaluations),  adaptive  class  structure,  and  the  flexibility  of  properties  to  be  of  any  class.  A  more 
thorough  description  of  these  capabilities  is  given  by  Veley.6 

The  GUI  that  was  developed  for  ASTROS  is  more  than  just  a  geometry  sketcher.  Properties  other  than  geometry  are 
defined  in  the  GUI  such  as  material  properties  and  loading  conditions  as  well  as  properties  such  as  design  variables  and 
design  constraints  that  are  used  to  define  an  optimization  problem.  Using  a  single  geometric  definition  and  the  non¬ 
geometric  information,  both  a  finite  element  model  and  an  aerodynamic  model  are  generated. 

At  the  current  stage  of  development,  the  GUI  does  not  define  a  complete  ASTROS  data  deck.  The  bulk  data  that  is 
missing  for  this  problem  are  those  entries  that  are  unique  to  the  modified  version  of  ASTROS  that  is  being  used  to 
handle  the  control  system  design.  Also,  a  comprehensive  generic  solution  control  object  which  the  user  can  tailor  to  suit 
the  user’s  model  is  still  under  development. 

Dependency  tracking  and  demand  driven  calculations  permit  parametric  relations  to  be  updated  automatically.  If  a 
change  is  made  to  the  sweep  angle,  span  or  chord  length,  the  ribs  and  spars  will  conform  to  the  new  planform. 
Similarly,  a  change  in  the  airfoil  shape  will  cause  the  height  of  the  ribs  and  spars  to  change  accordingly. 

The  adaptive  class  structure  of  AML  allows  for  ease  in  redefining  a  model.  Ribs,  spars  and  actuators  may  be  added  and 
deleted  from  the  model  with  a  few  clicks  of  the  mouse  button.  These  features  provide  much  ease  in  creating  new 
models  from  existing  models  and  significantly  reduces  the  time  to  develop  additional  models. 

4.  The  Flexible  Wing 

The  state-of-the-art  in  design  optimization  software  falls  somewhat  short  of  the  need  to  simultaneously  consider  all  the 
design  requirements  from  the  various  disciplines  connected  with  aircraft  design.  Thus,  obtaining  an  “optimum”  design 
in  an  integrated  fashion  is  not  yet  possible.  This  leads  to  the  design  process  consisting  of  a  combination  of 
optimizations  based  on  important  design  requirements  and  then  using  this  “optimized”  design  as  a  basis  for 
modifications  to  satisfy  other  design  or  performance  requirements  from  other  disciplines.  This  multiple  step  design 
approach  may  be  necessary  to  satisfy  global  and  local  requirements  form  different  disciplines. 

The  approach  used  in  this  effort  is  a  two  step  process.  First,  an  optimum  design  satisfying  strength  and  frequency 
distribution  requirements  is  obtained.  This  optimization  problem  is  solved  using  ASTROS  version  20  for  a  design 
condition  of  9-g  symmetric  pull-up  maneuver  at  M=0.85  at  30  psi.  Second,  using  a  modified  version  of  ASTROS,  the 
control  system  is  generated  to  twist  and  camber  the  optimum  wing  to  achieve  the  target  flexible  roll  rate. 

One  of  the  main  elements  of  the  flexible  wing  concept  is  that  rolling  moment  control  is  achieved  by  controling  the  wing 
twist  and  camber.  The  elastic  wing  twist  and  camber  needed  for  control  are  achieved  through  linear  actuators  mounted 
diagonally  across  the  internal  structure.  These  actuators  are  modeled  as  rods  exerting  tensile/compressive  loads  to  twist 
and  camber  the  wing  to  achieve  specified  roll  rates.  This  permits  assessments  between  energy  requirements  for  twist 
and  camber  control  to  achieve  the  required  roll  performance  at  various  dynamic  pressures. 

The  equilibrium  equations  for  steady  roll  maneuver  can  be  written  as7 

[K]{x]  +  <?[r]r[A]{a}  +  [F]{u}  =  0  (1) 

where  [F]  is  the  structural  stiffness  matrix  of  the  finite  element  model,  { x }  is  the  vector  of  nodal  displacement 
vector,  q  is  the  dynamic  pressure,  [TJ  is  the  transformation  matrix  from  structural  degrees  of  freedom  to  the 
aerodynamic  degrees  of  freedom,  {a}  is  the  vector  of  angle  of  attack  at  the  aerodynamic  panels,  [F]  is  the  applied 
actuator  load  distribution  matrix  and  {m}  is  the  vector  of  actuator  stimuli.The  product  [F]{m}  is  the  vector  of 
control  forces  generated  at  the  structural  node  points  due  to  the  actuator  forces  {u}  or  due  to  the  voltages  applied  to 
the  solid  state  actuators.In  the  latter  case  the  elements  of  matrix  [F]  would  depend  on  the  number  of  stacks,  number  of 
cycles  and  the  properties  of  the  solid  state  actuators  in  addition  to  the  direction  cosines  associated  with  the  actuators. 

The  displacement  vector  {x}  can  be  defined  as  a  linear  combination  of  rigid  body  modes  and  vibration  modes  as 


(2) 


{*}  =  +  =  MW 

where  [V,.]  is  the  rigid  body  mode,  [X)/^]  is  the  specified  number  of  antisymmetric  low  frequency  modes,  {T|r}  is  the 
generalized  rigid  body  displacements  and  in  present  case  it  is  equal  to  the  roll  angle  (j).  {T|  }  in  Equation  (2)  is  the 
elastic  displacements.  The  subscript  r  and  e  are  used  to  indicate  rigid  body  and  elastic  vibration  modes. 

The  angle  of  attack  {a}  at  the  control  points  of  aerodynamic  panels  can  be  written  as 

{«}  =  i[T]{x}  (3) 

where  V  is  the  free  stream  velocity.  Using  Equation  (2)  the  angle  of  attack  can  be  written  as 


{a} 


-[71 

V 


m  p  +  u 


d_Ve 

L  dx  J 


{lU 


where  p  is  the  roll  rate  (p  ,  U  is  the  chord  wise  component  of  V  in  the  direction  x , 
mode  gradient  with  respect  to  the  x  coordinate.  Equation  (4)  can  be  written  as 
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is  the  matrix  of  flexible 
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where  ’P  and  ’P  are  the  rigid  and  vibration  modes  expressed  in  aerodynamic  degrees  of  freedom.  Using  Equation  (5) 
the  equilibrium  Equation  (1)  can  be  written  in  terms  of  generalized  coordinates  as 


[K]jJL\  +  q[A]\JLl  +  [F][u]  =  0  (6) 


where  [A"]  = 
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Solving  Equation  (6)  for  the  generalized  coordinates 
gives 


=  tm  +  ^AirVuM},  do) 


The  roll  rate  p  is  given  by 


(11) 


P  =  (12) 

and  the  generalized  coordinates  of  the  elastic  modes  are  give 
by 


{r|J  =  [C]{u}  (13) 


For  the  specified  distribution  of  actuator  forces  {  u  }  Equation  Figure  3 .  Aerodynamic  model  of  first  wing 

(13)  can  be  used  to  calculate  the  roll  rate  p.  However,  for 
achieving  the  target  roll  rate  p  at  specified  Mach  number  and 

specified  dynamic  pressure,  the  optimal  control  design  approach  was  used  to  determine  the  distribution  of  actuator 
forces  necessary  to  deform  the  wing. 


5.  Application 

The  Adaptive  Modeling  Language  (AML)  Graphical  User 
Interface  (GUI)  developed  for  ASTROS  was  used  to 
model  two  similar  generic  fighter  wings.  The  finite 
element  model  developed  for  the  first  wing  which  is  more 
typical  of  fighter  wings  is  shown  in  Fig.  2  and  the 
aerodynamic  model  of  the  same  wing  is  shown  in  Fig.  3. 

The  semi-span  (fuselage  to  tip)  of  the  first  wing  is  180  in 
and  the  root  chord  is  163.9  in.  In  order  to  demonstrate  the 
effect  of  flexibility  in  the  wing,  the  second  wing  (shown 
in  Fig.  4)  is  a  modified  version  of  the  first  wing  in  which 
the  semi-span  and  root-chord  were  changed  to  be  246  in 
and  128.75  in,  respectively.  These  values  were  chosen  to 
maintain  the  same  wing  area  for  both  wings. 

In  order  to  assess  the  performance  of  the  AML-GUI  as  a 
model  generator  and  modifier,  a  timing  summary  was 
kept  of  each  of  the  activities  that  constitute  model  Figure  4.  Finite  element  model  of  modified  wing 

development.  The  timing  summary  is  given  in  Table  1. 

Not  included  in  the  summary  is  the  time  spent  reviewing  and  checking  the  entered  data.  The  total  time  in  front  of  the 
computer  for  defining  the  first  wing  was  approximately  2  hours  and  for  the  second  wing  approximately  one  half  hour. 
These  times  represent  the  time  of  someone  who  is  familiar  with  the  program  and  has  all  of  the  data  ready  to  enter  into 
the  program.  Note  that  the  amount  of  data  needed  to  be  computed  ahead  of  time  is  also  reduced  in  that  much  of  the  data 
is  “  translated”  from  one  form  to  another  within  the  program. 

The  functions  timed  are  as  follows:  The  definition  of  the 
planform  which  includes  the  span,  root  chord,  tip  chord, 
airfoil  shape,  twist  camber,  etc.  The  substructure  step 
includes  laying  out  the  drafting  lines  that  help  in  the 
definition  of  the  substructure,  and  the  placement  of  the 
wing  box,  ribs,  spars,  actuators  and  the  wing  attachment. 

The  properties  step  is  assigning  the  initial  sizes  and  the 
materials  (which  come  from  a  material  catalogue  that  is 
used  in  several  AML  programs)  to  the  wing  components. 

The  extra  points  and  boundary  conditions  (BCs)  include 
the  definition  of  additional  ‘‘drafting”  points,  RBARs, 

MPCs,  SPCs  and  SUPORT  conditions.  The  time  to 
generate  the  extra  points  is  approximately  the  same  as  the 
time  to  generate  the  boundary  conditions.  The  masses  step 
is  the  addition  of  non- structural  masses  to  the  model.  The 
aero  grid  step  is  the  definition  of  the  divisions  of  the 
planform  used  to  generate  the  aerodynamic  grid.  The  FE 
mesh  step  includes  the  time  to  generate  a  finite  element 
mesh,  define  design  constraints,  define  the  eigenvalue 
solution  method  and  write  out  the  ASTROS  data  deck. 

The  finite  element  mesh  is  completely  regenerated  each  time  and  is  for  the  most  part  a  push  the  button  and  get  a  cup  of 
coffee  type  function,  whereas  the  rest  of  the  development  time  is  nearly  constantly  hands  on  inputting  of  data. 

It  is  worthy  of  note  front  the  development  time  summary  that  the  primary  change  in  the  model  occurs  at  the  planform 
step  (a  conceptual  level  definition)  and  that  change  would  generally  affect  everything  that  follows.  However,  with  the 
dependent  relationships  between  each  function,  as  the  span  and  the  root  chord  change,  the  substructure,  masses  and 
aerodynamic  grid  move  with  it  and  do  not  need  to  be  redefined.  The  initial  properties  for  the  optimization  are  the  same 
and,  thus  do  not  need  to  be  changed.  Unfortunately,  the  current  method  used  to  define  the  extra  points  is  not  linked  to 
the  geometry  and  so  those  points  do  not  move  with  the  geometry  as  the  wing  is  stretched.  Consequently,  the  extra 
points  and  associated  boundary  conditions  need  to  be  redefined.  However,  those  boundary  conditions  that  are  not 
associated  with  the  extra  points  that  were  relocated  do  not  need  to  be  redefined.  Although  the  parameters  for  defining 
the  finite  element  mesh  are  dependent  on  other  properties,  the  mesh  itself  is  not  and  needs  to  be  regenerated  each  time. 
While  this  may  be  somewhat  time  consuming  for  large  meshes,  it  is  not  labor  intensive  so  the  engineer  may  perform 
other  tasks  while  the  finite  element  mesh  is  being  computed.  Although  no  relationship  changes  were  made  to  the 
substructure,  properties,  masses  or  aero  grid  in  this  example,  the  capability  does  exist. 


Table  1:  Timing  Summary  of  Model  Creation  (min:sec) 


Step 

Wing  1 

Wing  2 

Planform 

0:42 

0:12 

Substructure 

14:23 

— 

Properties 

3:33 

— 

Extra  Points/BCs 

24:50 

9:57 

Masses 

2:30 

— 

Aero  Grid 

1:07 

— 

FE  Mesh 

10:00 

10:00 

TOTAL 

58:00 

20:10 

The  two  wings  were  initially  designed 
for  a  9-g  symmetric  pull-up  maneuver 
at  Mach  0.85,  and  30  psi.  The 
optimized  designs  were  then  sent  to  a 
modified  version  of  ASTROS  with 
anti-symmetric  boundary  conditions 
to  generate  data  required  for 
designing  the  control  system  to 
achieve  a  120  degrees/sec  roll  rate. 

The  deformations  of  each  of  the 
wings  required  to  achieve  the 
specified  roll  rate  are  shown  in  Fig.  5 
and  Fig.  6. 

This  control  system  was  designed  for 
several  Mach  numbers  and  the 
corresponding  dynamic  pressure  for 
roll  maneuvers  at  20,000  ft.  The  strain 

energy  required  to  warp  the  wing  is  also  computed  and  is  given  in  Table  2. 

For  both  wings,  as  the  dynamic 
pressure  increases,  the  amount  of 
energy  required  to  warp  the  wing 
decreases.  This  is  as  expected  and 
demonstrates  how  the  flexibility  of 
the  wing  can  be  used  to  an  advantage 
for  maneuverability. 

Since  the  first  ten  frequencies  of  the 
wings  are  similar  (a  constraint  was 
placed  on  the  first  frequency  to  be 
greater  than  or  equal  to  3.0),  it  is 
concluded  that  the  stiffnesses  are  also 
similar.  The  energy  required  to  flex 
the  higher  aspect  ratio  wing  (wing  2) 
is  significantly  less  at  the  lower 
dynamic  pressures,  and  at  high 
dynamic  pressures,  the  lower  aspect 
ratio  wing  requires  less  energy  to  achieve  the  specified  roll  rate. 

It  is  believed  that  this  occurs  because  of  the  difference  in  roll  damping  between  the  two  wings.  Roll  damping  is  the 
effect  of  the  apparent  angle  of  attack  that  the  wing  sees  due  to  the  rolling  motion.  For  the  same  roll  rate,  flight  condition 
and  wing  area,  a  higher  aspect  ratio  wing  has  a  higher  roll  damping.  The  beneficial  effects  of  the  flexibility  of  the  high 
aspect  ratio  wing  are  competing  against  its  higher  roll  damping.  At  lower  speeds,  the  flexibility  benefits  outweigh  the 
roll  damping. 

6.  Summary 

The  linking  of  data  between  various  levels  of  analysis  using  AML  (e.g.,  linking  conceptual  design  parameters  to  the 
preliminary  design  models)  allows  data  decks  of  more  advanced  analysis  packages  to  be  generated  more  quickly  thus 
allowing  quicker  turn  around  on  the  higher  level  analyses.  If  the  time  to  modify  existing  models  is  significantly 
reduced,  the  effects  of  conceptual  level  decisions  on  more  detailed  designs  will  be  significantly  better  understood. 

While  this  effort  demonstrates  the  potential  for  quick  turn  around  of  finite  element  and  aerodynamic  models  for 
changes  in  planform  geometry,  the  AML-ASTROS  GUI  still  requires  much  work  to  enhance  the  robustness  of  the 
program  and  the  flexibility  of  the  program.  For  example,  the  program  is  currently  limited  in  the  planform  and  airfoil 
definitions  available.  The  definition  of  the  extra  points  and  boundary  conditions  could  be  modified  to  enhance  the 
dependent  relationships  so  that  they  do  not  need  to  be  redefined  when  the  planform  geometry  changes. 

This  paper  shows  the  benefit  of  being  able  to  quickly  incorporate  a  planform  change  into  the  finite  element  model.  The 
optimum  planform  and  substructure  design  for  roll  control  with  this  concept  are  linked.  The  substructure  design  effects 
the  stiffness  and  the  planform  effects  the  roll  damping.  The  design  process  to  exploit  this  concept  must  allow  for 
planform  and  substructural  design  simultaneously. 


Figure  5.  Deformation  of  Wing  1  needed  to  achieve  roll  rate 
of  120  degrees/sec  at  Mach  2.0 


Table  2:  Strain  Energy  Requirements  to  Achieve  120  deg/sec  Flexible  Roll  Rate  (Altitude 

20,000  ft) 


Mach 

Pressure 

(lb/in2) 

Strain  Energy a 
(ft  lbs) 

Strain  Energy15 
(ft  lbs) 

1.20 

6.81 

2427 

1537 

1.40 

9.27 

1945 

1173 

1.60 

12.1 

1519 

985 

1.80 

15.3 

1186 

862 

2.00 

18.9 

935 

775 

2.20 

22.9 

749 

709 

2.40 

27.2 

606 

657 

2.60 

32.0 

499 

580 

2.80 

37.1 

417 

560 

3.00 

42.6 

353 

550 

a.  Wing  1 

b.  Wing  2 
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